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Abstract 

The symmetry approach to the determination of Jacobi's last multiplier is inverted to provide a 
source of additional symmetries for the Euler-Poinsot system. These additional symmetries are 
nonlocal. They provide the symmetries for the representation of the complete symmetry group 
of the system. 
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1 Jacobi's Last Multiplier 

The method of Jacobi's last multiplier [[TJ, [TJ, [TJ (see also |§ [p 320, p 335 and pp 342-347] for 
a summary of these three papers of Jacobi) provides a means to determine an integrating factor, 
M, of the partial differential equation 



i=l 



(1) 



or its equivalent associated Lagrange's system 

dx\ dx2 



dx r , 



ai a 2 a n 

Provided sufficient information about the system ([!]) / (^) is known, the multiplier is given by 

d(f,u 1 ,u 2 , ■ ■ ■ ,uo n -i) 



(2) 



d(x 1 ,x 2 , ... ,x n ) 



MAf, 



(3) 



where 



d(f,ui,u 2 , • • • ,^ n -i) 
d(xi,x 2 , ... ,x n ) 
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and uji, . . . , io n -i are n — 1 solutions of (|l]) or, equivalently, first integrals of (|2|). As a consequence, 
one can prove that each multiplier is a solution of the partial differential equation 



i=l 



d(Ma t ) 



0. 



(5) 



A different combination of the integrals can produce a different multiplier, M'. The ratio M'/M 
is a solution of (|lj) or a first integral of (|2|) , which may be trivial as in the application of the Poisson- 
Jacobi Theorem |2l|, [ 10 1 to the determination of additional first integrals. A scholarly essay on 
the history of the Poisson-Jacobi theorem which Jacobi considered [1C] 



la plus profonde decouverte de M. PoissonF] 

and the pervasive influence of Jacobi's work upon Lie can be found in Q. 

In its original form the method of Jacobi's last multiplier required almost complete knowledge 
of the system under consideration. Since the existence of a first integral is consequent upon the 
existence of symmetry, one is not surprised that Lie [(j^] [pp 333-347] provided a symmetric route 
to the determination of Jacobi's last multiplier. A more transparent treatment is given by Bianchi 
|] [pp 456-464]. 

Suppose that we know n — 1 symmetries of (§|)/(0) 



l,n - 1. 



(6) 



Then Jacobi's last multiplier is also given by 



in the case that where now 



A = det 



ai, 

6i, 



M = - 



Ct2, 
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n-1,2, 



(7) 
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Jacobi's last multiplier provides an incestuous interrelationship between symmetries, first inte- 
grals and integrating factors for well-endowed systems. The practicality of this interrelationship was 
somewhat diminished in the past due to the effort required to evaluate the determinants of matrices 
of even moderate size. For example a Newtonian problem in three dimensions would require the 
evaluation of the determinant of a 6 x 6 matrix (7 x 7 if one considers time dependence, ie a% = 1). 
This possibly explains the omission of discussion of the method by postclassical authors such as 
Cohen 0, Dickson || and Eisenhart Q. The ready availability of computer algebra systems has 
rendered this method an attractive alternative, say, for the determination of first integrals given 
symmetries. To take a trivial example the 'free particle' has the Newtonian equation 



V = O 
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ll the most profound discovery by Mr Poisson'; see also the letter of Jacobi jnj 
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with the Lie point symmetries 



r 2 = td Ul + d U2 
T 3 = uid Ul + u 2 d U2 

r 4 = d t 

We find, for example, the determinants 



T 5 = td t - u 2 d U2 

T 6 = t 2 d t + tuid ui + (ui - tu 2 ) d U2 
T 7 = uid t - u\d U2 

T 8 = tu\d t + u\d Ul + (ui - tu 2 ) u 2 d, 
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(ui - tu 2 y 



(ui - tu 2 ) ul, 



tu\ u\ (u\ — tu 2 ) u 2 
using T 2 and T3; T3 and Tq; and T4 and Tg respectively, and the integrals 

A 2 

h = — = tu 2 - Ul 
^1 

r A3 2 

/2= aT = U2 



(10) 
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as expected. (Note that there are twenty-eight possible determinants to be calculated. Some of 
these are zero.) 



2 Complete Symmetry Groups 



In 1994 Krause |T6| introduced a new concept of the complete symmetry group of a system by defin- 
ing it as the group represented by the set of symmetries required to specify the system completely. 
There is not necessarily any relationship between the symmetries required to specify completely 
a system and its point symmetries. Thus Andriopoulos et al JTJ reported the complete symmetry 
group of the 'free particle' to be (D (B S T 2 ), the semidirect sum of dilations and translations in 
the plane, with the symmetries being three of the usual eight Lie point symmetries of (|| a), and that 
of the Ermakov-Pinney equation |27j to be ^3^ (sl(2,R)) with the symmetries being the three 
Lie point symmetries of that equation. On the other hand Krause fl~6|l reported that an additional 
three nonlocal symmetries are necessary to specify the Kepler Problem completely since the five 
Lie point symmetries of the three-dimensional Kepler Problem are insufficient to the purpose. This 
contains the implication that eight Lie symmetries are necessary to specify completely the Kepler 



Problem. However, a more careful analysis of sufficiency by Nucci et al [26|, based on the method 
of reduction of order proposed by Nucci in 1996 (2*|] and her interactive code for the determina- 
tion of Lie symmetries pl| , |2^ ], has revealed that the equation for the three-dimensional Kepler 
Problem is completely specified by six Lie symmetries with the algebra A\ © {A\ © s {2Ai © 2Ai}} 
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Ai(B{D © s {T2 © T2}}), where the subalgebra D(B S T2 is that associated with the first integrals 
of the one-dimensional simple harmonic oscillator |l9| (equally any second order ordinary differen- 
tial equation with the algebra sl(3,R)) to which the Kepler Problem reduces naturally under 
the method of reduction of order. 

One assumes that, when he devised his method of the last multiplier, the original intention of 
Jacobi was to determine integrating factors and that the adaptation from integral to symmetries 
by Lie was of like intention. However, (|8|) in combination with (||) suggests the possibility to 
determine symmetries provided the multiplier is known. The general solution of (||) is equivalent 
to the solution of @ unless one has the opportunity to perceive a particular solution without real 
effort. A particular case in point is when the functions di(x) are independent of x% for then (||) has 
the solution that M is a constant, taken to be chosen as a convenient value, which in this instance is 
not a 'trivial' solution. Since one now has an M, one can attempt to determine a further symmetry 
by solving (||) with one row of the matrix the coefficient functions of the unknown symmetry. One 
may infer that Lie and Bianchi had in mind point and contact symmetries in their treatments of 
Jacobi's last multiplier from the basis of symmetries rather than the original approach through 
first integrals used by Jacobi. However, as is common with many of the theoretical properties 
and applications of symmetries, there is no statement of the variable dependence of the coefficient 
functions in the method of Jacobi's last multiplier required for the method to hold. Consequently 
the considerations above apply equally to determination of nonlocal symmetries, in particular 
nonlocal symmetries of the type used by Krause, in which the nonlocality is found in the coefficient 
function of the independent variable, for autonomous systems. For such systems one of the known 
symmetries is dt which is represented in the matrix of (|8|) by the row (1,0,... ,0). In the Laplace 
expansion of the determinant the only possible nonzero terms must contain the first element of this 
row. If the unknown symmetry is 

r n = vd t + Gid Ui , (12) 

V does not appear in the expression for A, only the G{. These may be selected at will to satisfy 
the requirement that A = 1, or other suitable constant, and for each selection a V be computed 
through the invariance of the system of first order ordinary differential equations under the action 
of Tn\ the first extension of T n . In principle this would permit n symmetries to be determined. 
However, that presumes the independence of the Gi, i = l,n. This need not be the case if the Ui 
come from the reduction of an nth order scalar ordinary differential equation and the imposition 
of point symmetries is made at the level of the nth order equation. 



3 The complete symmetry group of the Euler-Poinsot system 

As an illustration of the ideas contained in §§1 and 2 we consider the simplest case of the motion 
of a rigid body which is the system governed by the Euler-Poinsot equations [||, 28] 

ui = — - — = W\ 
C — A 

io 2 = — — — w 3 wi = W 2 (13) 
A-B 



in which u := (u>i,uj 2 ,uj 3 ) t is the angular velocity and A, B and C are the principal moments of 
inertia. It is a commonplace that the system (|l~3| ) possesses respectively the two first integrals and 
Lie point symmetries 



E = \ [Aul + Buj 2 2 + Cujj) 
Ti = dt 



I? = A 2 co 2 + B 2 to 2 + C 2 u 2 
T 2 = -td t + uiidu. . 



(14) 



We note that Nucci's method of reduction of order [23] looks for first integrals in which one variable 
is missing J20l and provides p4j 



h = uj\B(A -B)- oj 2 3 C{C - A) 
I 2 = uj 2 C(B - C) - ujfA(A - B) 
h = uj 2 A(C -A)- uj 2 B{B - C) 



(15) 



as three conserved quantities, which are obviously not independent and can be constructed from 
the two independent integrals, E and L 2 , by the respective elimination of uj\, lu 2 and uj^ from them. 

We note that U{ is absent from the right hand side of d>i in (^) and so one solution of d|]) 
for Jacobi's last multiplier is a constant. The condition that T n ( |i~2| ) be a third symmetry of the 
Euler-Poinsot system (13) is 



A = Det 



1 Wi w 2 w 3 

10 

— t UJ\ u 2 uj 3 

V G\ G2 G3 



(16) 



where a is the constant to which we may assign some convenient value. We obtain 

AujihGi + Buj 2 I 2 G 2 + Cuj 3 I 3 G 3 = aABC. 

We take 

aBC 



Go 



0, G 3 







et eye 



(17) 



(18) 



and set a = I\/BC et eye to obtain the three sets of coefficient functions for the dependent variables 



— ,0,0 



0,— ,0 

ljj 2 



1 

0,0, — 

L0 3 



(19) 



We note that the sets in ( |l^ ) provide a basis and that other combinations could be taken. We keep 
the forms ( |l9|) simply for their present simplicity and find their subsequent utility. 

It remains to determine V. Consider the first set in (|l9|). The corresponding symmetry is 
written as 



r 3 = v 3 d t + — d Wl 

—2 + V3 ) - U2V 3 d^ 2 



iu 3 V 3 d L 



(20) 
(21) 
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and the action of on the Euler-Poinsot system, (|13"D, is 



lo 2 
"I 



C-A u 3 
B ~i 
A- Bu> 2 



(22) 



C ux 

in which (22b, c) are consistent with ( |i~3| ) and (|22|a) gives 

dt 



2 ■ 



A similar calculation applies to second and third of (p!9|). 
We obtain the three nonlocal symmetries 

dt 



(23) 



dt 
dt 



dt H 5 Wl 



a* + —a, 

U2 



UI2 



(24) 



d t + —5, 
w 3 



^'3 



for the Euler-Poinsot system (13) 



Since the system (|l^) is of the first order and autonomous, it is in a suitable form for reduction 
of order. We set y — U3 as the new independent variable. This is an arbitrary choice. Equally u>i 
or u>2 could be chosen as the independent variable. The results do not differ. The reduced system 
is 



dwi 
dy 

duj 2 



(B - C)C y 
A(A - B) ui 
(C - A)C y 



and inherits the symmetries 



dy B(A - B) oj 2 



f 2 = vid ui + uj 2 du 2 + yd y 



(25) 



r 3 = 
f 4 = 

?5 



-a 



1 
1 



-8y 

y 



(26) 



with the Lie Brackets 



r 2 ,r 3 



= -2r 3 , 
f 2 ,f 4 



r 3 ,r4 

= -2f 4 , 
f 2 ,f 5 



= 0, 

f 3 ,f 5 
2f 5 . 



r 4 ,r 5 
= 0, 



(27) 



G 



Consider a general two-dimensional system 

do>i 



dy 
dy 



f 2 {u 1 ,u 2 ,y), 



(28) 



of which the system (p5|) is a specific instance. We determine which of the four symmetries 
12,... ,T$ are necessary to specify (p5|) given (p8|). (This does beg the question of the appro- 
priateness of these four symmetries, but this will eventually become apparent.) 
The actions of 



fjf = LOid^ + u 2 d W2 + yd y 



= — d Ul 








— Uuj 2 


2°^9 


L0 2 


u4 2 



(29) 



on fl2q a) give 



dfi , dh dh 

y^r + w i« — I" — 

/i 1 d/i 



1 9/i 



(30) 



from which it is evident that 



U>2 0u>2 



ft(ui,u 2 ,y) = K 



y_ 

Ml 



(31) 



and hence (^ a) is recovered. 

The same symmetries acting on (p8|b) lead to (pB]b) and so the three symmetries, T2, T3 and 
T4 are a representation of the complete symmetry group of the system (f25|). By means of a similar 
calculation we see that the triplets of T2, T3 and T$ and of T2, and are also representations 



of the complete symmetry group of (25). 

The listing of Lie Brackets in (|27j) shows that the Lie algebra in each of the three cases is 
A\ (B s 2Ai or D © s T 2 , a representation of the pseudo-Euclidean group E(l,l) of dilations and 
translations in the plane. 



This is not be end of the story. Consider about the actions of T 3 
obtain three constraints on fi, videlicet 

h 1 dh 



r' ±J and f' 1] on (1281 a). We 



.,,2 



l_dh 

V dy 



(32) 
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when (]28| a) is taken into account. It is obvious that (^5| a) is recovered. Similarly (pqb) reduces 
to (|fb). 

Consequently the three symmetries T3, and T5 provide a representation of the complete 
symmetry group of (]25|). In this case the algebra of the symmetries is the abelian 3A^\ and not 
the A\ ® s 2Ai of the previous algebras. 

The system (25) is not the system (13) and one cannot expect that the complete symmetry 



group of (13) would be that of ( |25D although from the point of view of differential equations both 
are equally described in terms of a three-dimensional phase space. We consider the two triplets, 
videlicet T3, T4 and T§ and I?2, T3 and T4 (equivalently I?2, T3 and and I?2, and T5 as we saw 
above) and their respective actions on the general system. 

^1 = fl(t, UJl,U}2,U} S ) 

^2 = fz(t, UJl,UJ 2 ,UJ3) (33) 
^3 = f3{t,Ux,UJ2,UJ 3 ) 



to which class the Euler-Poinsot system ( |l3D belongs. 

The actions of the first extensions of I^, and videlicet 



if = ( - f - 2 ) d t + 1^, + ^ + (34) 



on (|33j a) are 



0J2 J ^2 W| u 2 

rl 1 ' = f — / — n J dt H <9(j 3 + —4^1 + —l^, 

\ ' w|y w 3 w| w| 



™ 3/l , 1 



h rp 9fi 1 dfi 

UJ2 Ot UJ 2 OLU2 

fi T dfi , 1 
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in which we have written T, = f dt/cof. The system ( |35| ) does not contain sufficient information to 
reduce (^ a) to the first of system (|l3|). The same applies for (j33|b) and (^3|c). The abelian group 
3Ai represented by r%, and 1^ is not the complete symmetry group of (|Tg|). 
The first extensions of T2 is 

T [ 2 1] = -td t + uid^ + u 2 duj2 + w 3^ 3 + 2cJ 1 <9< i , 1 + 2w 2 <^ 3 . ( 36 ) 



2 At first glance the existence of 3Ai for a second order system may seem to be at odds with the general result that 
a second order ordinary differential equation - a second order system - cannot admit "&A\ . However, ( |25[ ) cannot be 
written as a scalar second order ordinary differential equation. This emphasises the general point that a higher-order 
scalar equation may be reduced to a system of first order ordinary differential equations, but the reverse process is 
not always possible. 
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The actions of this, and rjj 4 on (|33|a) are 

2/i: 



dfx ^ dh ^ dh dh 
ac atJ! 0LV2 



From the first of @) 



nn 9fl 1 9/l 



fx = t 2 F 1 (u,v,w), 



(37) 



(38) 



where u, v and w are the three characteristics independent of fx, videlicet tux, tu2 and ^3. The 
substitution of p^) into (B%) and (57b) gives, respectively, 



1 dFi _ Ti 
u du t 3 
1 dFi Fx 

v dv v 2 



dFi dF 1 dFi 
-2F 1 + u^ 1 + v^- + w 



du 



dv 



dw 



T2 
t 3 



dF x dF x dF x 
-2Fx + u^- + v^- + w 



du 



dv 



dw 



(39) 



and now the situation is entirely different since the t dependence outside of the characteristics is 
isolated in the coefficients of the terms within crochets in both ( |39| a) and (|39|b). Consequently we 
have the three terms separately zero, ie 







1 dFx 

u du 

I^El-El = o 

v dv v 2 

dFx dFx dFx 
u— h v— h w- 



(40) 



du 



dv 



dw 



2Fx. 



We recover (|l3|a). Like calculations recover (||b) and (g|c). 

The complete symmetry group of the Euler-Poinsot system, @, is E(l, 1) (<^> D® S T 2 ). There 
are three equivalent representations. 



4 Discussion 

In this paper we have brought together several disparate ideas to arrive at something of a question 
mark. The symmetry-based version of Jacobi's last multiplier has been turned on its head, as 
it were, to provide a means to calculate nonlocal symmetries, in this instance, for the Euler- 
Poinsot system given the last multiplier. With these nonlocal symmetries we can identify the set of 
symmetries which specify completely the Euler-Poinsot system from all possible classes of systems 
of three first order equations. The number of symmetries required for this specification is three. 

This is surprising in terms of general expectations. The general first order equation for a system 
with one independent and n dependent variables is 

d>i = fi(t,ux, ■ ■ ■ ,u n ). (41) 
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With the application of each symmetry the number of variables in fi is reduced by one so that, 
after the application of n symmetries, there remains a general function of one characteristic. The 
application of a further symmetry specifies that general function. Thus we would expect that the 
complete symmetry group would have a representation in terms of an algebra of n + 1 elements. In 
fact one would regard an algebra of fewer than n + 1 elements as being quite unusual since it would 
imply an unexpected degree of connectedness amongst the coefficient functions of the different 
variables in the symmetries. (Such an example is found most dramatically in the case of the Kepler 
Problem @.) 



It is already known that a given system may possess more than one representation of its 
complete symmetry group. Indeed this is quite standard for any equation of the second order 
possessing eight symmetries. 

In the case of the Euler-Poinsot system we have found a rather intriguing result. In the reduced 
two-dimensional system, (|25|), we found that they were two inequivalent representations of the 
complete symmetry group, videlicet A\ ® s 2^4. 1 and "&A\. This lack of uniqueness does not persist 
when one returns to the three-dimensional system, properly known as the Euler-Poinsot system, 
for then at the latter symmetry group, videlicet 3^4i, falls away as a representation of the complete 
symmetry group. We do find that in conjunction with Y\ the three symmetries of ?>A\ do specify 
(|l~3|), but the four-dimensional algebra is not a candidate as a representation of the complete 
symmetry group since its dimensionality is not minimal. 

It would be interesting to find other examples of systems exhibiting similar properties^]. Cer- 
tainly the present result does place something of a question mark against the interpretation of the 
concept of a complete symmetry group as being the group of the symmetries which completely 
specify the equation although it does this without detracting from the inherent interest of the con- 
cept of a complete symmetry group. In fact one must seriously consider the identification of the 
characteristic system for a given problem, in this case whether it is the three-dimensional system 
fli~3|), which is the standard Euler-Poinsot system, or the reduced two-dimensional system (p5f). 
Curiously enough a similar question has arisen in the case of the Painleve Property for systems of 
first-order ordinary differential equations (Ti 
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